The line parameters and ratios as the physical probe of the line emitting regions in AGN 



D. Ilic a , J. Kovacevic b , L. C. Popovic b 

" Department of Astronomy, Faculty of Mathematics, University of Belgrade, Studentski trg 16, 1 1000 Belgrade, Serbia 
b Astronomical Observatory, Volgina 7, 1 1060 Belgrade, Serbia 



. . .Abstract 

on • 

f*^ Here we discuss the physical conditions in the emission line regions (ELR) of active galactic nuclei (AGN), with the special 
emphasize on the unresolved problems, e.g. the stratification of the Broad Line Region (BLR) or the failure of the photoionization 
to explain the strong observed optical Fe II emission. We use here different line fluxes in order to probe the properties of the ELR, 
(3J[)such as the hydrogen Balmer lines (Ha to He), the helium lines from two subsequent ionization levels (He II /14686 and He I /15876) 
and the strongest Fe II lines in the wavelength interval /1/14400 - 5400 A. We found that the hydrogen Balmer and helium lines can 
be used for the estimates of the physical parameters of the BLR, and we show that the Fe II emission is mostly emitted from an 
intermediate line region (ILR), that is located further away from the central continuum source than the BLR. 
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1. Introduction 



O r Active galactic nuclei (AGN) are a ubiquitous phenomena, 
' in a sense that most galaxies experience some sort of activity 
^ in their nucleus during their evolution. The most accepted sce- 
nario of the AGN structure is that they are powered by the ac- 
C"j cretion of matter onto a supermassive black hole (SMBH). One 
of the ways to study the inner emitting regions of an AGN, is by 
__ analyzing its emission lines, i.e. the broad (BELs) and narrow 
emission lines (NELs). So far many papers and textbooks are 
devoted to the physical p roperties of the emission line r egions 



(N 



(ELR) (see e.g. [ Boroson & Green, 1992; iSulentic et all 12000 



Osterbrock & Ferland, 2006, and references therein), but, there 



are still many open issues. 

Spectroscopy, in general, offers different methods for di 



agnostics of the em i tting plasma (see e.g. Griemi 1997 



f^**) [Osterbrock & Ferland , 2006 ). but these methods could not be 



properly used to probe the physical conditions in some ELR of 
• i-H AGN, as e.g. in the Broad Line Region (BLR), since the for- 
bidden lines, that are usually employed in plasma diagnostics of 
H the Narrow Line Region (NLR) or HII regions, are not present 



in the BELs spectrum. Particularly, it is difficult to find a direct 
method which would only use the observed BELs to determine 
the temperature and density in the BLR. On the other hand, the 
optical Fe II (AA 4400-5400 A) lines are one of the most inter- 
esting features in the AGN spectrum. The origin of the optical 
Fe II extreme emission and the geometrical place of the Fe II 
emission region in AGN, are still open questions. Also, there 
are many correlations of the Fe II lines and other AGN spectral 
properties which need physical explanation such as: EW Fe II 

vs - tmSt 'i ew Fe 11 vs - ro I 111 ' EW Fe 11 and FWHM 

of HjS, etc. dBoroson & GreenLll992h . 



In order to estimate the physical conditions (such as the tem- 
perature and hydrogen density) of the BLR we use the Balmer 
and helium line ratios obtained in two ways: (i) using the pho- 
toionization code CLOUDY, a spectral synthesis code designed 
to simulate conditions within a plasma and model the result- 
ing spectrum, and (ii) extracting a sample of AGN from the 
Sloan Digital Sky Survey (SDSS) database. We investigate 
these line ratios in order to find conditions in the BLR where so- 
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called Boltzmann-plot (BP) method is applicable (Griem, 1997; 
iPopovid, 120031 l2006al) . For these special cases, we study the 
correlations between the average temperature, hydrogen den- 
sity and He II/He I line ratio. Moreover, we present an inves- 
tigation of the optical Fe II emission in AGN, for which we 
have used an additional sample of 1 1 1 AGN from the SDSS 
database. The strongest Fe II lines are identified and classified 
into four groups according to the lower level of the transition: 
4 F, 6 S, 4 G and 2 D1. In this progress report, we report our re- 
cent investigations of the physical and kinematical properties 
of the BLR and Fe II emitting region. This report is organized 
as follows: in §2 we describe the numerical simulations of the 
BLR and briefly introduce the BP method, and give the anal- 
ysis of the simulated BELs; in §3 we study the SDSS sample 
of hydrogen Balmer and helium lines, while in §4 the selection 
and analysis of the SDSS sample of Fe II lines are given; in §5 
we discuss some results and finally, in §6 our conclusions are 
given. 

2. The BEL simulations 

In order to study the BELs, we simulated the BLR emis- 
sion line spectrum from different grids of the BLR photoion- 
ization mode l s using the CLOUDY code (version C07.02.01: 
Ferland et all 1 19981 iFerlandl [2006). Input parameters for the 
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simul ations are chosen to match the standard cond itions in the 
BLR dFerlancj 12005 iKorista & Goadll2000ll2004 . i.e. the so- 
lar chemical abundances, the constant hydrogen density, the 
code's AGN template for the incident continuum shape. We 
compute an emission-line spectrum for the coordinate pair of 
hydrogen gas density «H[crrT 3 ] and hydrogen-ionizing photon 
flux <l>H[cm _2 s _1 ]. The grid dimensions spanned 4 orders of 
magnitude in each direction, and with an origin of log «h = 8, 
log Oh = 17 was stepped in 0.2 dex increments. A column 
density jViitcni -2 ] was kept constant in producing one grid of 
simulations. Even though many authors claim that the most 
probable value of the BLR column density i s Nh = 10 23 cm~ 2 
dPumont. Collin-Souffrin & Nazarovai 1 19981; iKorista & Goad , 
l2000l 12004). we produce here 5 different grids of models chang- 
ing the column density in the range A^h = 10 21 - 10 25 cm~ 2 . 

We further analyze the BEL fluxes^ from the CLOUDY 
grids of models. We consider in our analysis the hydro- 
gen Balmer lines (Ha to He) and the flux ratio R of the he- 
lium lines He II ,14686 and He I ,15876, defined as R = 
F(HeIL14686)/F(HeL15876), where F is the line flux. We par- 
ticularly consider these helium lines, since they are the lines 
of the same element but in two different ionization states, thus 
their ratio He II ,14686/He I /15876 is sensitive to the change in 
the temperature and density (see e.g. Grieml 1997b . Besides, 
these lines are in the same spectral range as Balmer lines. Ad- 
ditionally, from the grids of models we consider in our analysis 
an averaged temperature, which is the electron temperature av- 
eraged over the BLR radius (T av further in the text). 



2.1. The Boltzmann-plot method for excitation temperature di- 
agnostics 

For the plasma of the length I that emits along the line of 
sight, assuming that the temperature and emitter density does 
not vary too much, the flux F u \ of a transition f rom the upper 
to the lower level (u -» D can be calculate d as (lGriemUl997 : 



Popovic, 2003, 2006a; P opovic et all 12008) 



he No 
Fui = —g u A u i— exp(-EJkT)( 

A L 

where A is the transition wavelength, g u is the statistical weight 
of the upper level, A u \ is transition probability, A^o is the av- 
eraged total number density of radiating species which effec- 
tively contribute to the line flux (which are not absorbed), Z is 
the partition function, E„ is the energy of the upper level, T 
is the averaged excitation temperature, and h, c and k are the 
Planck constant, the speed of light, and the Boltzmann con- 
stant, respectively. The additional assumption made here is 
that the population of the upper level in the transition follows 
the Saha-Bolt zmann distribution (for more detailed derivation 
check Popovic et all 2008b . From the above equation comes 
the so-called Boltzmann plot (BP) that can be used to estimate 
the excitation temperature r exc in the BLR 
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Figure 1 : Some examples of the BP applied on the Balmer line ratios calculated 
with the CLOUDY models. The normalized intensities F„ are calculated using 
the Balmer lines normalized on the H/3 flux. In the right corner of every plot, 
the BP temperature Tbp, the error of the BP fit / and the average temperature 
Tav are given. 



i°g W (F„) = log! 



Fui ■ A 

g U A„l 



= B-AE U 



The CLOUDY code gives all line fluxes normalized to the H/J flux. Since 
it has no influence in our analysis, we have used these values. 
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Figure 2: The error of the BP fit / in the hydrogen-density vs. ionizing flux plane for column densities Nu £ [10 21 - 10 25 ] cm 2 with contours inside which the BP 
fit error is smaller than 10%, 20% and 30%. The asterisks represent the simulations which satisfied three defined constraints (see text for details). 



where B and A are BP parameters, and A = log w (e)/kT exc « 
5040/ r exc is the temperature indicator. Therefore, for one line 
series (as e.g. Balmer line series) if the population of the up- 
per energy states (n > 3^3 can be described with the Saha- 
Boltzmann distribution, then by applying the last equation to 
that line series and obtaining the value of the parameter A, we 
can estimate the excitation temperature of the region where 
these lines are originating. We should emphasize that the ad- 
ditional assumption in the BP method is that the Balmer lines 
are originating in the same emitting region. 

In spite of the fact that the BP method has some obvious ad- 
vantages, i.e. it is only using the measured Balmer line fluxes, 
that are easily observed, to estimate the excitation temperature, 
one should take into account some possible problems that could 
occur when using the emission line to determine the BLR physi- 
cal properties: (i) the line profiles are usually very complex and 
indicate that more than one component contribute to the total 
line flux, therefore one should be aware of the multi-component 
BLR structure when estimating the BEL fluxes; (ii) the Balmer 
lines do not have to necessarily originate in the same region, 
e.g. there some indications that the Ha and HB line are forming 
in two distinct regions since it ha s been shown that the H/3 l ine is 
systematically broader than Ha (IShapovalova et all 12008): (iii) 
there are different mechanisms contribution to the line forma- 
tion. The photoionization seems to be working well, but other 



2 We note here that since the emission deexcitation goes as u — » / it is not 
necessary that level / has the Saha-Boltzmann distribution. 



heating mechanisms should be taken into account. 

2.2. The analysis of the simulated BELs 

The first step in analyzing the BEL ratios is to apply the BP 
method on the Balmer line ratios given by CLOUDY to esti- 
mate the parameter A, from which we then calculate the BP 
temperature, i.e. the excitation temperature (7bp further in the 
text) of the region where Balmer lines are formed. Also, from 
the best-fitting of the normalized line ratios, we obtain the error 
of the BP fit (/ further in the text). A few examples of the BP 
applied on the Balmer lines simulated with CLOUDY for col- 
umn density of Nn = 10 23 cirT 2 are presented in Fig.Q] In many 
cases a satisfactory BP fit is not obtained, i.e. / has pretty large 
values. This is more noticeable in the case of higher values of 
the hydrogen density and ionizing-photon flux, hence we plot 
the error of the BP fit / in the hydrogen density vs. ionizing 
flux plane for all 5 grids of models of different column density 
in Fig. |2] (only the contours inside which / is less than 10%, 
20% and 30% are given). From Fig.|2]can be seen that if the BP 
method could be considered valid if / is less than 10% (even- 
tually 20% in the measured spectra) the parameter space where 
this is valid is pretty constrained for all column densities A^h- 

Out of all photoionization models, we select just those that 
follow these constraints: (i) the error of the BP fit / less than 
10%, (ii) the average temperature r av less than 20,000 K, as for 
larger temperatures the Balmer line ratios are not sensitive to 
the temperature changes and the BP method cannot be applied 
Popovic 2006b), and (iii) the ratio of helium lines R less than 



Table I: The best-fitting parameters A,- and B, of the function T av = A,- + B, ■ R, 
and Cj and D, of the function log rin = D, / (C, +R) for different column densities 





logA^H 


A; 


Bi 


Di 
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[cm -2 ] 


[K] 


[K] 


[cm -2 ] 




21 


4422±689 


10468+1697 


5.04+0.22 


0.97+0.05 


22 


3308+572 


10413±1098 


3.76+0.17 


0.70+0.04 


23 


3486+200 


7116+288 


2.91+0.10 


0.51+0.03 


24 


4634±87 


3326+134 


2.37+0.10 


0.40±0.03 


25 


5208+119 


1899+152 


2.16+0.12 


0.41+0.04 



2. Following this criteria, we constructed 5 samples of simu- 
lated BELs, of single column density, labeled as e.g. CD21 for 
column density Nh = 10 21 cirT 2 , etc. The models that satisfy 
these 3 constraints are represented with asterisks in Figure [2] 



2.3. The correlations between the BEL ratios and the BLR 
physical properties 

We investigate in more details the five sets of simulated spec- 
tra defined above. First, we plot the average temperature r av of 
the emitting region, one of the outputs of the model, with re- 
spect to the ratio R of the helium lines for all 5 sets (Fig. [3] 
upper). We fit data sets with the linear function T av = A + B ■ R 
and the best fitting results are given in Table Q] Another possi- 
ble connection of the BLR physical parameters and BELs, can 
be obtained from the relation between the hydrogen density and 
the helium line ratio R. We plot the logarithm of the hydrogen 
density and the ratio R (Fig. [3] bottom), and we fitted it with 
the function: \ognn - D/(C + R), where is given in the 
units of 10 7 cm~ 3 . It can be seen from Fig. [3] (bottom) that even 
when the column density changes, the same dependance of the 
hydrogen density on the helium line ratio R remains. The best 
fitting results are given in Table Q] 

3. The hydrogen Balmer and helium line sample 

We compare our results of the numerical simulations with 
the observed data. For that we consider our measurements of 
well defi ned sample of 90 AG N taken from the SDSS spectral 



database (La Mura et al., 2007), for which the Blamer line ra- 



tios have been precisely measured and the temperature param- 
eter A has been estimated using the BP method applied on the 
Balmer line series (see for details La Mura et all 120071) . We 



found that for this sample in approx 50% of cases, the BP 
method can be applied. In this analysis we consider objects 
that have error of the BP fit smaller than 2O9S0 and the BP tem- 
perature smaller than 20 000 K. This reduced our sample to 48 
objects. 

For that sample, we measure the fluxes of the helium lines He 
II /14686 and He I /15876, being particularly careful in the case 
of He II /14686 line, where there is contamination from the Fe II 
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3 Even though in numerical simulations the error was taken to be smaller 
than 10%, here in real measurements, we take larger range of fitting errors due 
to the additional observational and line flux measurement errors. 



0.5 1 1.5 2 

He II (X4686) / He I (15876) 

Figure 3: The average temperature T dv of the emitting medium (upper) and the 
hydrogen density (bottom) as the function of the ratio R of the He II /14686 
and He I A5816 lines, for five selected sets of simulated spectra of different 
column density Nh = 10 2 ' - 10 25 cirT 3 . The dashed line represents the best- 
fitting with the function T- dv = A, + B, ■ R for the average temperature and 
log «h = DjliCj + R) for the hydrogen density. The name of the set is specified 
on each plot. 
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4475 4550 4625 4700 4775 4850 4925 5000 5075 
Wavelength (in A) 

Figure 4: An example of the Gaussian decomposition of the Hy8 region in the 
case of the object SDSSJ1300+5641. The Gaussian components are shown 
below the observed spectrum, where the helium line He II .14686 is fitted with 
two Gaussians. The residual spectrum is given in the bottom. The Fe II template 
is presented with dashed line. 



multiplet and nearby H/3 line. The local continuum have been 
subtracted from the spectra, as well as the cont ribution of the 
host galaxy (see for details lLa Mura et al .1 120071 for details). 
For the He II /14686 line measure ments, we perform Ga ussian 



decomposition (for details see e.g. iPopovic et aU 120041) in or- 



der to subtract the contribution of the Fe II multiplet, H/3 and 
narrow He II line (as an example see Fig. @J. A standard de- 
viation is taken as an error of the line flux measurements. The 
cases when helium lines could not be measured are not consid- 
ered (e.g. when lines were too noisy or the contribution of Fe II 
was too strong and could not be properly subtracted), thus our 
sample is reduced to 20 objects. 

Using the above relations between the physical properties of 
the BLR and helium line ratio R (Table [TJ, we estimate the 
average temperature and the hydrogen density of the BLR of 
these AGN assuming different column densities. We obtain 
the following ranges of physical parameters in the BLR: av- 
erage temperature T- dv = 5500 - 17600 K and hydrogen density 
n H = 10 8 ' 3 - 10 IL6 cirT 3 . 

Moreover, in Fig. [5] we plot the ratio of the helium lines as a 
function of the BP temperature T^p for the SDSS sample. As 
can be seen in Fig. [5] there is a weak correlation (correlation 
coefficient r = 0.50, po - 0.02321) between BP temperature 
and the He line ratio. This is also in agreement with the BLR 
physics, where for higher temperatures one should expect to 
have stronger He II than He I lines. We should note here that if 
we exclude the point that is clearly high above the best-fitting 
line (see Fig. |5j, the correlation coefficient is slightly better r = 
0.62 (po = 0.00473). 



4. The Fe II line emitting region 

We selected a sample of 1 1 1 AGN from Sloan Digital Sky 
Survey (SDSS) according to the following criteria: high signal 
to noise ratio (S /N > 20), good quality of the pixels, negligible 
contribution of the stellar component and a good coverage (near 
uniform) of redshifts from to 0.7. 
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Figure 5: The He II lines ratio as a function of the BP temperature 7bp for the 
SDSS sample. The linear best-fitting is presented with solid line. 



As result, our sample contained 111 spectra, from which 58 
have all Balmer lines, and the rest 53 are without the Ha line 
(due to cosmological redshifts). After removing the influence 
of the Galactic extinction and underlying continuum, the Fe II 
lines in the region of the H/3 line are fitted with the calculated 
Fe II template. The template consist of the 33 strongest Fe II 
lines in the AA 4400-5400 A range, separated into four groups 
according to the lower level of transition: 4 F, 6 S, 4 G and 2 D1. 
Since all considered Fe II lines probably originate in the same 
region, values of their shift and width are assumed to be the 
same in one AGN, while the intensities are different. The [O 
III] /U4959, 5007 A, He II ^4686 A, [N II] AA6548, 6583 A 
and Balmer lines (H/3 and Ha) are fitted with a sum of Gaus- 
sians. We assumed that Balmer lines have three components: 
the narrow, intermediate and very broad, that correspond to the 
narrow (NLR), intermediate (ILR) and very broad (VBLR) line 
region, respectively. The NLR components in one AGN are as- 
sumed to have the same width and shift as narrow [O III] and [N 
II] lines. The examples of the fit of the H/3 and Ha line region 
are shown in Fig. [6] and Fig. [7] respectively. 




Fe II (U 4472-4731 A) 



Re II (14924 A) : 

[ Fell (1501 8 A) Fe II (U 51 69-5325 A) 



4475 4550 4625 4700 4775 4850 4925 5000 5075 5150 5225 5300 5375 
Wavelength (in A) 

Figure 6: The spectrum of SDSS J075101.42 + 29174000.00 in the AA 4400- 
5500 A range: A) observed spectra (dots) and the best fit (solid line). B) H/3 
is fitted with the sum of three Gaussians which represent emission from NLR, 
ILR and VBLR. The [O III] AA4959, 5007 A lines are fitted with two Gaussians 
for each line of the doublet and He II /14686 A is fitted with one broad Gaussian. 
Template of Fe II is denoted with a dashed line, and represented also separately 
in panel C (below). 
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6480 6500 6520 6540 6560 6580 6600 6620 6640 6660 
Wavelength (in A) 

Figure 7: Example of the fit of the SDSS J075 101 .42 + 29174000.00 spectrum 
in the AA 6400-6800 A range: Ha is fitted with the sum of three Gaussians 
which represent emission from NLR, ILR and VBLR and [N II] ,U6548, 6583 
A lines are fitted with one Gaussian for each line of doublet. Narrow emission 
lines (Ho- NLR and [N II]) are denoted with dashed lines. 

In order to investigate the geometrical place of the Fe II emis- 
sion region in AGN, we analyzed kinematical connection be- 
tween the Fe II and Balmer emission region by comparing the 
line shifts and widths, obtained from the best-fitting. We as- 
sume that the broadening of the lines is caused by the Doppler 
effect from the random velocities of emission clouds, while the 
shift of the line is caused by the systemic motion of the emis- 
sion gas. The results are presented in Fig. [8][9j 

The widths of the Fe II and H/3 (upper) and Ha (bottom) 
lines are compared in Fig. [8] The X-axis gives the Fe II line 
widths, while the Y-axis gives the parameters of the widths of 
the H/3 (upper) and Ha (bottom) components: the NLR (trian- 
gles), ILR (circles) and VBLR (squares). Dotted lines show the 
averaged values of widths: Fe II (vertical line) and H/3 and Ha 
line components (horizontal lines). On the sample of the 111 
AGN that contain the H/3 line (upper), the average value of Fe 
II widths is 1770 + 580 km/s, while the average values for H/3 
components are: 310+ 160 km/s (NLR), 1850 + 660 km/s (ILR) 
and 4690 ±1210 km/s (VBLR). Considering the selected sam- 
ple of 58 AGN which have the Ha line (bottom) the averaged 
value of the Fe II width is 1770 ± 570 km/s and for the Ha 
components: 270 + 130 km/s (NLR), 1470 ± 620 km/s (ILR) 
and 4070 ± 1080 km/s (VBLR). It is obvious that the averaged 
Fe II width (1770 km/s) is very close to the averaged widths of 
the ILR component of both the H/3 and Ha line (1850 km/s and 
1470 km/s, respectively), while the averaged widths of the NLR 
and VBLR components are significantly different. 

Fig. [9] shows the correlation between the Fe II widths and 
the H/3 ILR ones (r = 0.67, p Q < 0.0001). Considering the 
Ha ILR component, correlation is stronger, and it is r = 0.72, 
po < 0.0001. We have also found that a weaker correlation 
between widths of the Fe II and Ha VBLR component (r = 
0.55, po < 0.0001) is present. Relations between the shifts of 
Fe II and Ha and H/3 (NLR, ILR and VBLR components) are 
also investigated. It is noticed a positive correlation with the 
shift of the Ha ILR component (r = 0.51, p Q < 0.0001). There 
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Figure 8: The widths of the Fe II lines compared to the widths of Hf} com- 
ponents (upper, the sample of 111 AGN) and Ha (bottom, the sample of 58 
AGN); On X-axis are widths of Fe II, and on Y-axis are the widths of the NLR 
(triangles), ILR (circles) and VBLR (squares) components of H/J (Ha). Dotted, 
vertical line shows the average value of Fe II widths, while dotted, horizontal 
lines show the average values of the Hfi (Ha) components. The average value of 
Fe II lines is the same or very close to average width value of ILR components 
of H/3 and Ho-. 



is no a significant correlation between the Fe II shift and other 
H a and H/3 components. 

Boroson & Greenl (|l992|) found negative trend between the 



equivalent width (EW) of Fe II and [O III] lines (r = -0.39), as 
well as the anticorrelation between the EW Fe II and EW ([O 
III])/EW (H/3) (r = -0.53). Thus, we tested the same relations 
on our sample of 1 1 1 AGN. We found slightly lower correla- 
tions then the previous authors: for EW Fe II vs. EW [O III] 
the correlation coefficient is r — -0.26 (po = 0.007) and for EW 
Fe II vs. EW ([O III])/EW (H/3) it is r = -0.37 (p < 0.0001). 

Then, we divided our initial sample of 1 1 1 AGN (0<z<0.7) 
in 6 subsamples, which are made by gradually removing the 
objects with low redshift in 6 steps. Subsamples contain object 
with redshifts in intervals: 0.1<z<0.7, 0.2<z<0.7, 0.3<z<0.7, 
0.4<z<0.7, 0.5<z<0.7 and 0.6<z<0.7. Relation between EW 
Fe II and EW [O III] (also EW [O III]/EW H/3) are analyzed 
separately for each subsample (see Table |2). We found that 
negative trend (r = - 0.37, P<0.0001) observed between EW Fe 
II and EW [O III]/EW H/3 in initial sample (0<z<0.7), progress 
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Figure 9: The correlation between the Fe II widths and widths of the H/3 ILR components (left), Ha ILR components (middle) and Ha VBLR components (right). 
In all cases, correlations are observed (0.67, 0.72, 0.62), indicating kinematical connection between these emission regions. 



Table 2: Correlations between equivalent widths (EW) of the Fe II and [O III] 
([O III]/H/?) lines for subsamples within different redshift intervals. Relation 
is fitted with function Y = A + B ■ X. The coefficient of correlations r, the 
corresponding measure of the significance of correlations pa , as well as the A 
and B coefficients, are shown in the table. In the last column, the number of 
objects in each subsample is given. 



VS. 




EW [O III] 


EW [O III]/EW HB 


No 


EW Fell 


r 


-0.26 


-0.37 


109 


(0<z<0.7) 


Po 


0.007 


< 0.0001 




E W Fell 


r 


-0.30 


-0.40 


92 


(0.1<z<0.7) 


Po 


0.004 


< 0.0001 




E W Fell 


r 


-0.42 


-0.50 


74 


(0.2<z<0.7) 


Po 


2.2E-4 


< 0.0001 




EW Fell 


r 


-0.44 


-0.56 


58 


(0.3<z<0.7) 


Po 


6.0E-4 


< 0.0001 




EW Fell 


r 


-0.56 


-0.62 


43 


(0.4<z<0.7) 


Po 


< 0.0001 


< 0.0001 




EW Fell 


r 


-0.64 


-0.70 


28 


(0.5<z<0.7) 


Po 


2.37E-4 


< 0.0001 




EW Fell 


r 


-0.68 


-0.72 


14 


(0.6<z<0.7) 


Po 


0.007 


0.004 





to significant anticorrelation (r~ - 0.70, P<0.0001) as we are 
removing low-redshift objects. The redshift dependance of the 
EW Fe II vs. EW [O III] anticorrelation may indicate some 
cosmological repercussion of the Fe II emission in AGN or de- 
pendence on luminosity. 

5. Results and discussion 

In this progress report we give some results of our recent 
investigations of different broad line parameters (the flux ratios, 
EW, FWHM, etc.) in order to probe the physics of the ELR, 
i.e. the hydrogen Balmer lines (Ha to He), the helium lines 
from two subsequent ionization levels (He II /14686 and He I 
/15876) and the strongest Fe II lines in the wavelength interval 
^4400 - 5400A. 

From the analysis of the hydrogen Balmer lines, we can say 
that the BP method gives valid results if the error of the BP fit 
/ is less than 10% (eventually 20% in the measured spectra). 



The parameter space of CLOUDY simulations where / < 10% 
is well constrained and in a similar range for different column 
densities A^h (Fig. This area of parameters contains lower 
ionizing fluxes and higher hydrogen densities, but depending 
on the column density the area increases keeping the same trend 
between the nu and <$>n, Therefore, for those parameters the BP 
method could be applied for the estimation of the excitation 
temperature. This indicates that for some physical conditions, 
even if we have the photoionization as the main heating pro- 
cess in the BLR, the hydrogen Balmer lines are produce in such 
way that they obey the Saha-Boltzmann equation, i.e. the BLR 
plasma is at least in the Partial Local Thermodynamical Equi- 
librium. 

The physical parameters of the simulations for which the 
BP method could be applied (/ < 10%) follow some rela- 
tions even when the column density changes. There is a lin- 
ear relation between the average temperature T av of the emit- 
ting region and helium lines ratio R: T m — A + B ■ R (Fig. |U 
Table [TJ. Also, the hydrogen density depends on the helium 
line ratio R as: log«// = D/(C + R), where tin is given in the 
units of 10 7 ctrT 3 (Figure |H Table Q]). For lower column den- 
sities (Nh = 10 21 - 10 22 cm~ 2 ) this linear trend for J av seems 
to be broken after some value of R (Fig. HI upper), therefore 
the results obtained for this column densities should be taken 
with caution. The ranges of the average temperature and hy- 
drogen density, obtained when these relations are applied to the 
observed sample of AGN taken from the SDSS database, are 
in good agreement with the previous estimates of the physical 
conditions in the BLR (lOsterbrock & Ferlandl |2006j). Having 
in mind problems given in §2.1 (i.e. different emitting regions 
of helium and Balmer lines or the multicomponent origin of 
BELs) the relations given in Table Q] could be use as a rough 
estimate of the BLR physical parameters from direct measure- 
ments. 

On the other hand, from the analysis of the Fe II emission, 
we found the significant correlation between the Fe II and Ha, 
HB ILR widths which indicate the ILR origin of the Fe II emis- 
sion. Also, the Fe II emission region is mainly characterized 
with random velocity of ~1800 km/s, that corresponds to the 
ILR origin. This result is in favor to the results obtained by 
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Pop ovic et al. ( 20041) . and with the recent findings that the op- 
tical Fe II line forming region seems to not be the same as th e 



BLR dKuehn et alll2008l:lHu et all koOS: Pop ovic et al.[|2009h 



Moreover, we found that the degree of anticorrelation of EW 
Fe II - EW [O III] (EW [O III]/H0) is redshift dependent, i.e. 
anticorrelation increases by selection of the sample with high 
redshift objects. 

6. Conclusions 

We used here the emission lines to estimate the physical pa- 
rameters of the BLR and properties of the Fe II line forming 
region. From our investigations we can outline the following 
conclusions: 

(i) In the case of the pure photoionized plasma, there can 
exist the conditions that the excitation of the hydrogen 
Balmer lines is sensitive to the temperature. In this case, 
the He II/He I line ratio can be used for the estimates of the 
averaged temperature and hydrogen density. This can be 
very useful for determining the physical properties of the 
BLR, and for a small sample of AGN we found that the 
BLR temperatures are ranging from ~ 5000 K to 18000 
K, and hydrogen densities from 10 s cirT 3 to 10 12 cm -3 . 
These values of the t emperature and density are expected 
in the BLR plasma dOsterbrock &Ferlandil2006l) . 

(ii) The optical Fe II emission is likely coming from a region 
that corresponds to the ILR. Also, there is a weak correla- 
tion with the VBLR component of Balmer lines indicating 
that a fraction emitted in the Fe II lines can originate also 
in the BLR. On the other hand, the correlation between 
the EW Fell and EW [ODI] (EW rOIIIl/EW Hy3) discov- 
ered bv lBoroson & Green (1992) seems to be sensitive on 
redshift (or luminosity) of AGN. 

Finally, the BLR physics (and geometry) still remains an 
open question, and we hope that with this work we approached 
some conclusions which will help us in understanding the BLR 
physics. 
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